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IMPROVING  JET  REACTOR  CONFIGURATION  FOR  PRODUCTION  OF  CARBON 

NANOTUBES 

ALEX  POVITSKY* 


Abstract.  The  jet  mixing  reactor  has  been  proposed  for  the  industrial  production  of  fullerene  carbon 
nanotubes.  Here  we  study  the  flowfield  of  this  reactor  using  the  SIMPLER  algorithm.  Hot  peripheral  jets  are 
used  to  enhance  heating  of  the  central  jet  by  mixing  with  the  ambiance  of  reactor.  Numerous  configurations 
of  peripheral  jets  with  various  number  of  jets,  distance  between  nozzles,  angles  between  the  central  jet  and 
a  peripheral  jet,  and  twisted  configuration  of  nozzles  are  considered.  Unlike  the  previous  studies  of  jet 
mixing,  the  optimal  configuration  of  peripheral  jets  produces  strong  non- uniformity  of  the  central  jet  in  a 
cross-section.  The  geometrical  shape  of  reactor  is  designed  to  obtain  a  uniform  temperature  of  a  catalyst. 

Key  words,  jet  mixing,  SIMPLER,  jet  heating,  turbulent  jets,  nanotube  production,  fullerenes 

Subject  classification.  Fluid  Mechanics 

1.  Introduction.  The  C^o  molecules  combined  in  nanotubes  have  unusual  and  very  useful  electronic 
and  mechanical  properties  [7],  [9].  However,  the  mass  production  of  nanotubes  still  presents  a  challenge. 
This  report  describes  the  CFD  modeling  and  optimization  of  a  jet  mixing  reactor  for  such  a  production.  The 
single-walled  carbon  nanotubes  may  be  produced  from  CO  with  the  help  of  a  catalyst.  Catalyst  particles 
get  into  the  reaction  chamber  via  the  central  CO  jet.  The  particles  must  be  heated  from  200 C  (central 
jet  temperature)  to  temperatures  close  to  1000C  (reactor  temperature)  as  rapidly  as  possible,  since  slowly 
heated  catalytic  particles  quickly  loose  their  catalytic  ability. 

The  mixing  reactor  studied  here  consists  of  an  arrangement  of  jets  designed  to  heat-up  catalytic  particles 
carried  by  a  turbulent  CO  jet  issuing  from  a  round  nozzle.  This  jet  entrains  the  hot  gas  within  the  reactor 
cavity  and  heats  the  jet  material.  However,  this  heating  is  considered  too  slow  to  produce  a  significant  amount 
of  carbon  nanotubes.  To  increase  mixing,  hot  turbulent  jets  issuing  from  several  nozzles  surrounding  the 
central  nozzle  are  used.  These  jets  are  inclined  to  the  central  jet  to  further  increase  the  mixing  rate.  The 
temperature  of  the  hot  jets  is  equal  to  the  average  reactor  temperature. 

The  distance  between  the  central  and  the  peripheral  nozzles,  the  number  of  peripheral  jets,  and  the 
angle  between  the  central  jet  and  a  peripheral  jet  are  taken  as  optimization  parameters.  Figure  1  shows  a 
typical  jet  interaction  geometry. 

The  behavior  of  a  single  jet  in  a  static  ambiance,  the  influence  of  the  density  difference  between  the  jet 
and  surrounding  gas,  the  effects  of  co-flowing  round  jets,  jet  rotation,  and  the  development  of  jets  in  a  cross- 
flow  stream  have  been  extensively  studied  [1],  [2]  and  results  are  summarized  as  empirical  formulas,  confirmed 
by  experiments  and  computations.  However,  the  three-dimensional  mixing  of  jets  at  high  intersection  angle 
a  includes  mutual  dependence  of  the  above  factors.  Therefore,  a  detailed  CFD  investigation  and  discussion 
of  the  physics  of  interacting  jets  is  needed  for  the  optimization  of  the  reactor. 

Holdeman  et  al.  [4]  studied  mixing  of  a  set  of  side  jets  with  a  flow  in  a  cylindrical  duct  and  obtained 
the  optimal  number  of  jets  for  a  given  jet-to-mainstream  momentum  flux  ratio  (or  mass  flux  ratio).  In  this 
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work,  they  consider  a  confined  and  uniform  central  flow  in  a  duct,  whereas  in  our  reactor  a  cold  central  jet 
is  injected  into  a  hot  ambient  gas.  In  the  current  study,  the  size  of  the  reactor  cross-section  is  larger  than 
that  of  the  central  nozzle,  therefore,  the  central  jet  is  more  similar  to  a  free  jet  than  to  a  duct  flow.  Thus, 
a  peripheral  jet  should  not  be  mixed  with  the  reactor  hot  gas.  Instead,  the  peripheral  jet  must  reach  the 
central  jet  keeping  its  axial  velocity  as  large  as  possible  and  then  cause  the  rapid  mixing  of  the  central  jet 
with  the  ambiance.  To  achieve  this  goal,  the  strongest  possible  peripheral  jets  are  used.  If  the  total  mass 
flux  of  peripheral  jets  is  fixed,  we  show  that  the  minimum  number  of  peripheral  jets  ensures  fast  mixing  of 
the  central  jet  with  the  reactor  gas.  Our  calculations  show  that  two  peripheral  jets  in  the  ( x ,  z)  plane  cause 
the  central  jet  to  spread  out  in  the  (y,z)  plane  leading  to  the  fast  heating  of  the  cold  central  jet  by  mixing 
with  hot  reactor  gas. 

To  further  increase  the  mixing  of  the  central  jet,  the  peripheral  jet  nozzles  were  drilled  off  the  z  axis  to 
cause  rotation  of  the  central  jet.  Aerodynamic  twisting  has  been  exploited  in  energy  technology  [6]  to  increase 
turbulent  transport.  However  when  heating  is  to  take  place  at  short  axial  distances,  a  nozzle  arrangement 
with  non-zero  total  angular  momentum  of  peripheral  jets  is  more  advantageous  to  the  basic  configuration 
due  to  the  diagonal  elongation  of  the  central  jet  cross-section  than  due  to  its  rotation. 

For  cases  where  a  is  greater  than  90°,  the  nozzles  are  installed  on  the  side  walls  of  the  reactor. 

The  requirements  of  higher  a,  short  distance  between  a  peripheral  jet  nozzle  and  its  intersection  with 
the  central  jet,  and  have  enough  space  in  a  cross-section  to  spread  the  central  jet  lead  to  optimization  of 
reactor  shape  in  this  study.  In  our  future  work,  we  shall  consider  such  an  optimization  for  non- Cartesian 
geometry  of  reactor. 

The  paper  is  organized  as  follows.  In  Section  2,  we  conduct  a  study  for  the  optimal  angle  a,  based  on 
empirical  expressions  for  jets  in  cross-flow.  This  is  considered  as  a  preliminary  estimation  before  the  CFD 
analysis.  In  Section  3,  we  present  a  CFD  framework  for  our  numerical  analysis.  In  Section  4,  we  study 
the  jets’  interaction  numerically  in  order  to  define  the  optimal  values  of  afore-mentioned  parameters  of  the 
jet  reactor.  In  Section  5,  we  explore  a  jet  arrangement  with  non-zero  total  angular  momentum  leading  to 
aerodynamic  twisting  of  the  central  flow,  and  the  counter-flow  jets’  arrangement  where  the  nozzles  are  drilled 
through  the  side  walls. 

2.  Simplified  model  of  mixing  of  jets.  Downstream  of  the  nozzles,  the  injected  cold  and  hot  jets  are 
assumed  to  develop  as  free  jets  until  they  intersect.  We  assume  that  these  jets  intersect  at  a  single  point  at 
the  centerline.  After  the  intersection,  the  central  jet  axis  remains  straight  due  to  the  symmetric  arrangement 
of  the  peripheral  jets,  whereas  each  of  the  hot  jets  is  deflected  from  its  original  direction  by  the  cross  flow  of 
the  central  jet. 

The  goal  is  to  ensure  fast  heating  of  a  central  cold  jet  by  peripheral  jets.  The  external  part  of  the  central 
jet  has  been  heated  by  mixing  with  the  hot  ambiance  before  the  intersection  with  peripheral  jets.  The 
problem  is  to  heat  up  the  area  of  the  central  jet  near  the  centerline.  Therefore,  we  consider  a  penetration 
length  (i.e.,  axial  distance  between  beginning  of  the  jet  interaction  and  the  intersection  point  of  jets’  axis)  as  a 
reasonable  indicator  of  fast  heating.  We  investigate  the  penetration  length  as  function  of  the  jet  intersection 
angle  a. 

For  a  round  turbulent  jet,  the  centerline  velocity  distribution  is  given  by  [1] 

(1)  Um  =  10.7  u0  — ^ —  if  X  >  Xi 

X  +  Xo 

and  the  distribution  of  velocity  in  a  cross-section  is  as  follows 

(2)  ux  —  um{  1  —  3?7  +  2r/  ), 
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where  77  —  y/b,  b  —  (x  4-  xq )tanfl,  tanfi  —  0.22,  =  2. 067*0,2;*  =  S.95ro,  um  is  the  local  velocity  at  the  jet 

axis,  Xi  is  the  length  of  initial  part  of  the  jet,  u0  is  the  jet  velocity  (x  =  0),  r0  is  the  jet  radius,  and  0  is  the 
jet  expansion  angle. 

Eqs  (1)  and  (2)  are  used  to  compute  the  jet  profiles  prior  to  their  interaction. 

Demuren  [2]  uses  an  empirical  model  for  the  trajectory  of  single  turbulent  jet  injected  normally  into  a 
cross  flow 


(3) 


w^ere  5^7  >  are  coordinates  of  the  jet  centerline,  the  empirical  coefficients  are  chosen  as  a  =  0.85, 

b  =  0.47,  and  c  =  0.36 [2],  the  momentum  flux  ratio  q  —  index  00  denotes  cross  flow,  and  index  m 

denotes  parameters  of  the  impinging  jet  at  point  A. 

In  our  case  the  central  cold  jet  remains  straight  due  to  the  symmetric  orientation  of  heating  jets.  Thus, 
we  can  consider  the  cold  jet  as  a  cross  flow.  However,  the  central  jet  is  not  equivalent  to  a  uniform  cross 
flow.  Thus  the  linear  Taylor  expansion  of  the  above  expression  is  used  to  compute  the  coordinates  of  the 
axis  of  the  peripheral  jet  at  discrete  points  (y(i),ar(i))  : 


(4) 


%eff 


where  q  ^(^00(7),  Poo(O)- 

In  the  above  equation,  any  of  the  heating  jets  is  assumed  to  be  a  circular  turbulent  jet  injected  normally 
into  the  cross  flow.  Peripheral  jets  are  no  longer  uniform  while  they  hit  the  central  jet.  The  effective  diameter 
of  a  peripheral  jet  before  it  impinges  into  the  central  jet  is  not  defined.  A  reasonable  guess  is  deff  =  d0. 
That  suggests  consideration  of  a  stream- line  tube  with  the  above  mentioned  diameter.  To  compute  <7,  only 
the  y-  component  of  a  peripheral  jet  is  taken  into  account.  The  computations  of  jet  trajectory  start  at 
point  A  and  last  until  the  jet  intersection  at  point  I  (see  Figure  1). 

The  optimal  mixing  of  jets  is  assumed  to  be  the  case  where  the  peripheral  jet  reaches  the  centerline  of 
the  central  jet  as  fast  as  possible,  i.e.  deviation  from  its  initial  direction  is  small.  Clearly,  this  is  not  the  case 
for  small  a,  where  jets  are  almost  parallel.  Even  for  large  a,  the  dense  central  jet  deviates  light  peripheral 
jets. 


Figure  2  shows  the  penetration  length,  normalized  by  the  orifice  radius,  versus  a.  These  computations 
lead  to  a  —  90°  which  corresponds  to  the  minimum  penetration  length.  However,  one  should  take  into 
account  the  large  hydraulic  resistance  that  occurs  for  large  a  and  consider  the  slope  of  the  curve  in  Figure  2 
rather  than  its  minimum  value.  The  slope  is  large  for  a  <  50°  and  small  for  larger  values  of  a.  Thus,  taking 
a  «  50°  is  a  reasonable  compromise  between  speedy  heating  and  complications  of  large  values  of  a. 

This  section  offers  useful  but  preliminary  analysis  of  jets’  mixing.  The  CFD  results  (Sections  4  and 
5)  show  that  the  central  jet  develops  strong  angular  non-uniformity  in  the  (y,z)  plane.  Obviously,  this 
important  process  in  mixing  is  not  taken  into  account  here.  Also,  the  CFD  results  supply  us  with  a  flowfield 
temperature  distribution  which  is  a  direct  measure  of  the  quality  of  heating. 


3.  CFD  model. 


% 


3.1.  Governing  equations.  The  Navier-Stokes  equations  of  gas  dynamics  with  an  appropriate  model 
of  turbulence  describe  the  mixing  of  jets.  The  system  of  transport  equations  is  given  by 


(5) 


-E 


djpUjF) 

dxi 


+  Sf  ~  0, 
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where  ar1,x2,  and  x3  are  Cartesian  coordinates,  F  =  Ui,  C/2,  f73,  fc,  e,  C,T  are  the  main  variables,  Ui,U2,U3 
are  velocity  components,  k  and  e  are  kinetic  energy  of  the  turbulence  and  the  turbulent  dissipation,  C  is  the 
mass  concentration  of  material  in  the  central  (cold)  jet,  T  is  the  temperature,  SF  is  the  source  term,  and  T 
is  the  transport  coefficient.  The  Standard  k  —  e  model  [8]  is  used  in  this  study  to  predict  turbulent  transport 
when  the  turbulent  viscosity  and  transport  coefficient  are  given  by 

r  — 

PrF ’ 

(6)  fieff  =  C„pk2/€, 

where  cp  =  0.09  is  the  coefficient  of  the  k  —  e  model  of  turbulence,  and  PrF  is  the  Prandtl  number  for 
variable  F. 


Table  1 

Prandtl  numbers  and  source  terms  for  governing  equations 


Variable 

Prandtl  number 

Source  Term 

Ui 

1. 

QP 

dxi 

c 

.75 

- 

T 

.9 

- 

k 

1. 

Q-ep 

€ 

1.3 

(< Ci <5  -  c2ep)e/k 

In  Table  1,  C\  =  1.44  and  c2  =  1-92  are  coefficients  of  the  h  —  e  model  of  turbulence,  and  Q  is  the  rate 
of  turbulence  energy  generation: 


(7) 


Q  —  P'eff 


\fajJ  \dxj  dxi 


For  the  special  case  where  Mach  number  M  <  0.3,  the  density  is  assumed  to  be  independent  of  pressure.  For 
mixing  of  chemically  inert  jets,  Eq.  (5),  with  zero  source  term,  is  solved  for  the  concentration  C  of  material 
in  the  cold  jet.  To  satisfy  the  boundary  conditions,  the  concentration  C  is  set  equal  to  unity  at  the  central 
nozzle  and  equal  to  zero  at  the  peripheral  nozzles.  Using  the  known  concentration  field,  the  local  density  is 
computed  by 


(®)  P  —  C Pcold  (1  C)photy 

where  pcoia  and  ph0t  are  densities  of  the  central  and  peripheral  jets. 

The  catalyst  particles  considered  are  small  (less  than  20 p)  in  terms  of  their  size  and  mass,  therefore, 
they  are  assumed  to  move  with  zero  velocity  relative  to  the  gas.  Thus,  the  same  concentration  field  C  shows 
the  spatial  distribution  of  catalyst  particles. 

The  temperature  field  is  computed  by  Eq.  (5)  with  zero  source  term,  where  the  boundary  conditions  are 
similar  to  those  for  the  concentration  field.  Still,  the  coefficients  of  the  PDE  for  temperature  are  different 
from  those  for  concentration  as  their  Prandtl  numbers  are  different.  Note  that  the  temperature  field  does 
not  affect  other  governing  equations. 

3.2.  Numerical  algorithm.  We  discretize  the  governing  PDEs  (5)  using  a  finite  volume  method  and 
a  structured  numerical  grid.  The  current  version  of  the  developed  computer  code  includes  discretization  on 
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uniform  or  non-uniform  grids  in  Cartesian  coordinates.  For  simplicity,  we  describe  the  uniform  grid  version 
here.  The  basic  grid  is  used  for  discretization  of  T,  C,  P,  K ,  and  e  governing  variables: 

xi (i)  =  {i-  l)hu  i  = 

X2 O')  =  -L2  4  (j  -  l)h2,  j  =  1, N2j 

(9)  x3(k)  =  -L3  +  (k  -  l)/i3,  k  —  lj  ...,N3y 

where  hi  =  Li/(iVi  -  1),  h2  =  2L2/(N2  -  1),  /&3  —  2L3/(N3  -  1)  are  grid  steps,  and  Pi  x  2 L2  x  2L3  is 
the  computational  volume.  The  symmetry  of  solution  is  one  indicator  of  the  correctness  of  computer  code, 
therefore,  we  do  not  halve  the  computational  domain  by  imposing  the  symmetric  boundary  conditions  unless 
it  is  stated  otherwise.  Staggered  grids  are  used  for  discretization  of  velocity  components: 

xu(i)  =  0.5(#i (i)  4-  Xi(i  —  1)),  i  —  2, ..., Ni 
xv(j)  =  0.5{x2(j)  +  x2(j  ”  1)),  j  =  2, ..., N2 

(10)  xw(k)  =  0.5(x3(fc)  +x3(k-  1)),  k  =  2,  ...,7V3, 

where  xw(l)  =  Xi(l)  -  0.5/ii,xi>(l)  =  x2(l)  -  0.5h2,xw(l)  —  x3(l)  -  0.5/i3,  xu(Ni  H-  1)  =  x\(Ni)  + 
0.5hi,xv(N2  4-  1)  =  x2(iV2)  -f  0.5h2,xw(N3  4-  1)  =  x3(N3)  4-  0.5/i3.  Grids  {xu,x2ix3},{x i,xv,x3}7  and 
{xi,x2,xw}  are  used  for  discretization  of  ?7i,t/2,and  U3y  respectively.  Basic  and  staggered  grids  in  a  2- 
D  section  are  shown  in  Figure  3.  The  finite-volume  method  for  discretization  of  Eq.  (5)  on  any  of  the 
described  above  numerical  grids  leads  to  the  following  set  of  equations  with  respect  to  unknown  values  of  F 
in  neighboring  grid  points: 

(11)  apFp  —  ajvFjv  +  asF$  4  &eFe  4  awFw  4-  apFp  4  (ibFb  4  Sf, 
where 


ap  =  dE  4  O'W  +  aN  4”  a>s  +  &t  +  &b  4 

(  (pUt)n-(pU2).  ,  (pU»)t~(pUz)b. 

Xlse  xltw  X2,n  X2)S  X3)£  X3j& 

1  r 

aE  = - max(-(pUi)e ; - 0.5(pPi)e;  0); 

Xl,e  -  Xltw  xl,E  ~  Xi,p 

1  r 

(12)  aw  - - max((pUi)w ; - h  0.5(pE/i)w;0). 

xl,e  xl,w  X\ X\ }W 

Here  indices  P,  W,  iV,  5,  T,  P  correspond  to  the  immediate  neighbors  of  the  point  P  (see  Figure  4).  Velocities 
in  corresponding  directions  are  computed  at  mid-points  e,  w ,  n,  s,  t,  b.  Formulas  for  coefficients  a#,  aw ,  ay,  a# 
are  similar  to  above  formulas  for  coefficients  aw  and  05.  This  form  of  discrete  coefficients  implies  a  central 
difference  second-order  scheme  while  |Pe|  <  2  and  a  first  order  upwind  scheme  otherwise  [5],  [3],  where 
Pe  =  pUAx/T.  To  compute  coefficients  (12),  the  turbulent  viscosity  and  density  at  the  point  e  are  calculated 
as  the  arithmetic  average  of  its  values  at  points  E  and  P.  Eqs.  (11)  with  coefficients  similar  to  Eqs.  (12)  are 
used  to  compute  velocities  on  staggered  grids. 

For  computations  of  the  pressure  field  and  correction  of  velocities,  that  are  made  to  satisfy  the  continuity 
equation,  Patankar’s  SIMPLER  (Semi-Implicit  Method  for  Pressure-Linked  Equations  Revised)  is  used  [5]. 
The  reduced  formulas  (11)  for  velocities  lead  to  the  following  expressions  connecting  the  velocity  and  pressure 
corrections: 

(13)  U'e  =  -  PB~pP  f 

(xi,j E  —  Xl,p)oe 
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where  ae  is  the  coefficient  ap  for  velocity  U\.  The  SIMPLE-type  algorithms  substitute  the  above  formulas 
in  the  continuity  equation  to  form  Eq.  (11)  for  the  pressure  correction.  The  set  of  discretized  equations  for 
pressure  correction  is  given  by 

(14)  bpPp  =  b^PN  4-  bsPg  4*  b&PE  4-  bwPw  4*  bpPp  4-  +&b-Pb  4-  5p' » 

where 


(15) 


bp  =  6jv  4-  6s  4-  6#  4-  b\v  4-  bp  4-  bp) 
1  Pe 


6  at  = 
6s  = 


£l,e  -  XiiW  dN(Xl,E  -  X\ ,p)  ’ 

1  p(*  -  1/2) 

^l.e  —  %l,w  ds(xltp  -  Xitw)  ’ 


-( 


w.  -  w. 


*El.< 


‘  Xl,w 


,  Wn-m.  ,  (pDi,)t  -  (pDaJj 

4 - 1 - 


X2  e  X2 


%3,t  ~  X3 tb 


)■ 


After  computing  the  P'  field,  velocity  corrections  are  computed  by  Eq.  (13).  SIMPLER  employs  the 
pressure-correction  equation  only  for  the  task  of  correcting  the  velocities  and  not  for  the  computation  of  the 
pressure  field. 

To  compute  pressure,  Eq.  (14)  is  used  with  coefficients  similar  to  (15),  where  pseudo-velocities  are 
used  instead  of  velocities  [5].  At  each  external  iteration  by  all  governing  equations,  the  Poisson  Eqs.  (14) 
for  pressure  and  for  pressure  correction  should  be  solved  up  to  a  high  level  of  convergence  to  ensure  the 
satisfaction  of  the  continuity  equation. 


4.  Arrangement  of  nozzles  at  the  reactor’s  entrance.  The  reactor  is  considered  as  a  box  with  rect¬ 
angular  entrance  section,  equipped  with  nozzles,  impermeable  side  walls,  and  open  exit.  Physical  variables 
are  normalized  by  the  central  nozzle  radius,  the  central  jet  speed  and  its  density  at  the  nozzle. 

The  cold  central  jet  and  hot  peripheral  jets  have  initial  temperatures  Tjet  =  200 C  and  Tenv  =  1000 C, 
respectively.  In  all  cases  considered,  the  nozzle  gas  velocity  of  a  peripheral  jet  is  equal  to  that  of  the  central 
jet.  In  Case  7  (two  peripheral  parallel  jets),  the  radius  of  the  central  orifice  and  each  of  peripheral  orifices 
are  equal.  For  a  fair  comparison,  the  total  mass  flux  and  initial  velocity  of  the  peripheral  jets  are  taken  the 
same  for  all  considered  configurations  of  nozzles.  Therefore,  the  radius  of  a  peripheral  nozzle  is  computed 
by: 


(16)  rv~  rp0/(0.5Ar?cos(a))0'5, 

where  rv o  is  the  radius  of  a  peripheral  orifice  in  Case  7. 

Input  data  for  the  considered  arrangements  of  nozzles  are  gathered  in  Table  2.  In  fact,  the  difference 
between  cases  is  fully  determined  by  the  boundary  conditions  at  the  entry  section.  The  boundary  conditions 
at  the  entry  cross-section  and  a  section  of  the  numerical  grid  in  (y,z)  plane  are  shown  in  Figure  5. 

The  quality  of  heating  of  the  catalyst  particles  is  estimated  from  the  particle  temperature  distribution 
at  the  reactor  exit  cross-section.  This  distribution  shows  the  heating  rate  on  average  and  the  uniformity  of 
heating.  The  height  of  a  histogram  column  is  computed  by 

(17)  hi=  f  plhCdxdy 

Js 

where  S  is  the  subset  of  the  cross-section  such  as  0j_i  <6  <6i,  and  6  =  (T  -  Tjet)/(Tenv  —  Tjet). 
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Table  2 

Arrangements  of  nozzles 


Case 

Distance,  r/rc 

Angle,  a 

Number  of  jets,  Nj 

1 

- 

- 

single  jet 

2 

2.0 

5 

4 

3 

2.0 

60 

4 

4 

10.0 

60 

4 

5 

10.0 

30 

2 

6 

10.0 

60 

2 

7 

10.0 

0 

2 

The  first  case  corresponds  to  a  single  jet  heated  up  by  spreading  into  a  hot  ambiance.  This  case  shows 
the  basic  rate  of  heating  for  the  simplest  arrangement  of  the  reactor  (see  Figures  6  and  11a).  The  second 
and  the  third  cases  represent  the  compact  arrangement  of  four  peripheral  jets  taken  as  close  as  possible  to 
the  central  jet.  The  former  case  has  almost  a  parallel  arrangement  of  jets,  whereas  in  the  latter  case  the 
angle  between  the  central  and  a  peripheral  jets  is  equal  to  60°.  Some  computational  flow  fields  in  the  (x,z) 
plane  are  shown  in  Figure  7.  The  co-flowing  peripheral  and  central  jets  lead  to  elongation  of  the  high-speed 
centerline  core  of  the  central  jet  in  comparison  to  the  single  jet  (compare  Figures  7a  and  6).  This  leads  to 
slower  mixing  of  the  central  jet  with  the  ambiance.  Although  the  co-flowing  component  of  the  peripheral 
jet  speed  is  halved  in  the  latter  case,  the  jet  elongation  is  still  noticeable  in  this  case  (see  Figure  7d).  The 
concentration  and  temperature  fields  show  the  behavior  similar  to  that  for  velocity  (and  momentum)  fields 
(see  Figure  7b,c,e,f). 

In  turn,  this  leads  to  slow  mixing  of  the  central  jet  with  the  hot  gas  and  the  exit  temperature  distribution 
of  particles  is  either  lower  than  that  for  the  single  jet  (Case  2)  or  only  marginally  higher  (Case  3)  (compare 
Figure  11a,  b,  and  c). 

To  reduce  the  co-flowing  effect,  the  peripheral  jet  nozzles  were  placed  at  some  distance  from  the  central 
jet  (Case  4).  This  leads  to  a  definite  increase  of  the  final  particle  temperatures  (compare  Figure  lid  to  Figure 
11a).  However,  the  distance  between  the  peripheral  jet  nozzles  and  the  jets’  intersection  should  be  small  in 
order  to  keep  the  peripheral  jet  speed.  As  it  was  mentioned  in  the  previous  Section,  the  core  velocity  of  the 
free  jet  remains  equal  to  its  nozzle  velocity  while  the  axial  distance  from  the  nozzle  is  less  than  10ro  —  12ro 
and  then  it  decreases  proportional  to  this  distance. 

To  increase  the  mass  and  momentum  flux  of  a  single  peripheral  jet,  configurations  with  two  peripheral 
jets  are  considered  (Cases  7,  5,  6)  and  computational  fields  are  shown  in  corresponding  Figures  8,  9,  and  10. 
In  Case  7  (parallel  jets),  the  inclination  of  the  peripheral  jets  towards  the  central  jet  occurs  only  due  to  the 
entrainment  of  these  light  hot  jets  by  the  dense  central  cold  jet.  The  difference  between  the  jets’  plane  and 
the  middle  plane  (compare  Figure  8a, c  to  b,d)  is  slight  and  the  quality  of  heating  is  close  to  that  for  Case 
1.  Therefore,  the  entrainment  of  light  (hot)  jets  by  the  dense  (cold)  jet  should  not  be  considered  as  a  major 
factor  of  the  cold  jet  heating. 

The  angular  variation  of  field  isolines  in  (y,z)  plane  becomes  considerably  more  prominent  for  initially 
inclined  jets  (Cases  5  and  6).  The  central  jet  narrows  in  the  plane  (x,z),  where  the  peripheral  jets  hit  it, 
and  spreads  aside  in  the  plane  (#,?/)  (see  Figures  9  and  10).  As  opposed  to  the  previous  cases,  the  heated 
particles  occupy  the  entire  exit  cross-section  (see  Figures  9d  and  lOd).  The  quality  of  heating  is  best  in  Case 
6  and  second  best  in  Case  5  (see  Figure  lle,f). 
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5.  Offset  and  side- wall  jet  configurations.  In  this  section  we  consider  peripheral  jets  with  non¬ 
zero  total  angular  momentum  (Case  8),  and  peripheral  jets  with  a  >  90°  (Case  9)  (see  Table  3).  Two 
peripheral  jets  are  taken  here  according  to  results  shown  in  the  previous  Section.  In  order  to  produce 


Table  3 

Offset  and  side-wall  configurations  of  peripheral  nozzles 


Case 

Coordinates  of  nozzles  (x,  y>  z) 

Angle,  a 

Cross-section  (y,  z) 

8 

(0,±r0,±10r0) 

60* 

20ro  x  20ro 

9 

(17.3rO;0,±20ro) 

120 

20ro  x  20ro 

10 

(ly.Sro.O.ilOro) 

120 

20ro  x  10ro 

angular  momentum  in  the  central  flow,  the  peripheral  nozzles  are  drilled  at  distance  ±lro  from  the  z  axis. 
In  this  case,  the  angle  a  is  the  angle  between  a  peripheral  jet  and  a  horizontal  plane  that  includes  the 
central  cold  jet.  This  configuration  of  jets  leads  to  slightly  better  mixing  than  that  in  Case  6  (see  Figure 
12b).  The  reason  for  better  mixing  is  the  resulting  cross-sectional  elongation  of  the  central  jet  along  the 
diagonal  of  the  cross-section  (see  Figure  12a).  as  opposed  to  elongation  in  the  y  direction  as  in  Case  6.  The 
simple  geometric  observation  that  the  square  diagonal  is  \/2  times  larger  than  the  square  edge  explains  the 
difference  in  mixing  between  Cases  8  and  6. 

For  x  <  30ro,  the  rotation  of  the  central  jet  remains  small  and  cannot  lead  to  a  substantial  increase  of 
jet  mixing. 

Next,  we  investigate  Case  9,  where  the  jet  intersection  angle  a  =  120°  between  the  central  and  a 
peripheral  jets.  To  make  it  feasible,  the  peripheral  nozzles  are  drilled  through  the  side  wall.  For  fair 
comparison  to  Case  6,  the  jets  are  designed  to  meet  at  the  same  axial  location.  The  higher  angle  between 
jets  causes  better  mixing,  however,  this  nozzle  arrangement  increases  the  distance  a  peripheral  jet  travels 
before  intersecting  the  central  jet.  From  geometric  consideration,  this  distance  in  Case  9  is  twice  as  large 
as  that  in  Case  6  (compare  Figure  13a  to  Figure  10a).  The  tradeoff  between  higher  angle  of  incidence  and 
smaller  maximum  local  peripheral  jet  velocity  may  lead  to  different  conclusions  about  the  efficiency  of  such 
an  arrangement  of  jets.  Here,  the  temperature  distribution  of  particles  is  slightly  different  from  that  in  Case 
6  (compare  Figure  Ilf  to  Figure  13d).  In  Case  9,  the  amounts  of  well-heated  particles  (0.90  <  9  <  0.95)  and 
poorly  heated  particles  (0.7  <  6  <  0.75)  increase  at  the  expense  of  intermediately  heated  particles. 

To  achieve  a  smaller  distance  between  a  peripheral  jet  nozzle  and  the  intersection  point  (Case  10,  Table 
3),  the  cross-section  of  reactor  is  changed  from  a  square  to  a  rectangle  in  such  a  way  that  its  height  was 
halved  and  its  width  remains  the  same  as  in  Case  9.  Other  parameters  in  Case  10  remain  the  same  as  in 
Case  9.  This  cross-sectional  shape  combines  the  strength  of  peripheral  jets  (as  their  traveling  distance  is 
halved)  and  provides  the  available  space  for  spreading  the  central  jet  (as  the  2/— edge  remains  unchanged). 
For  these  cases  we  impose  symmetric  boundary  conditions  at  the  symmetry  planes  and  compute  a  quarter 
of  the  computational  volume  (see  Figure  14).  The  impact  of  the  peripheral  jet  on  the  central  jet  is  stronger 
than  that  in  Case  9  (compare  Figure  14a  to  Figure  13a).  Some  particles  mix  with  the  reactor  ambiance 
upstream  of  the  jets’  intersection  (Figure  14b)  because  the  strong  peripheral  jets  act  as  an  obstacle  to  the 
central  jet  and  create  stagnation-type  flow  at  the  centerline.  The  mixing  of  particles  within  the  reactor 
downstream  of  the  jets’  intersection  occurs  through  the  impact  of  the  crosswise  component  of  flow  on  the 
side  walls  and  the  creation  of  vortices  in  the  cross-section  shown  in  Figure  14c.  Intensive  mixing  leads  to  an 
even  heating  of  particles  (see  Figure  14d). 
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6.  Conclusion.  Fast  mixing  and  heating  of  the  central  jet  occurs  in  cases  with  strong  angular  non¬ 
uniformity  of  the  merged  jet  in  a  cross-section  (yyz)  and  wide-spreading  of  particles.  This  appears  to  be  a 
key  process  for  the  aerodynamic  mixing  of  the  central  jet  with  the  reactor  ambiance.  Unlike  previous  studies 
with  high  angular  uniformity  in  a  cross-section  and  use  of  multiple  peripheral  jets,  the  maximum  strength  of 
each  peripheral  jet  suggests  a  configuration  of  two  peripheral  jets  with  high  intersection  angle  between  the 
central  and  a  peripheral  jets.  The  optimal  distance  between  peripheral  nozzles  and  the  central  jet  results 
from  the  tradeoff  between  the  strength  of  a  peripheral  jet  while  it  hits  the  central  jet  and  the  co-flowing 
effect  delaying  mixing  of  the  central  jet. 

Arranging  peripheral  nozzles  in  such  a  way  that  the  total  angular  momentum  of  the  peripheral  jets 
is  non-zero  helps  to  control  the  direction  of  spread  of  the  central  jet  and  may  increase  its  mixing.  Using 
counter-flow  jets  increases  the  spread-up  of  the  central  jet  and  requires  detailed  investigation  of  the  reactor 
geometry. 

Further  studies  will  be  conducted  to  investigate  the  effect  of  cross-section  shape  and  size  on  mixing. 
Possible  improvements  include  variable  axial  size  of  cross-section.  A  small  effective  diameter  of  the  reactor 
cross-section,  where  peripheral  jets  hit  the  central  jet,  provides  the  strength  of  peripheral  jets,  whereas 
larger  reactor  cross-section  downstream  provides  the  necessary  space  for  elongation  of  the  central  jet’s  cross- 
section.  Perhaps,  an  elliptical  cross-section  provides  the  required  elongation  length  in  one  direction  and  a 
short  distance  between  the  peripheral  nozzles  and  the  central  jet  in  the  other  direction.  To  perform  this 
study,  the  wall  boundary  layer  model  and  the  ability  to  handle  complex  geometries  will  be  incorporated  in 
the  computer  code. 
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Fig.  3.  Staggered  grid  (2-D  ( x,y )  section) ,  O  15  numerical  grid  for  T,C,P,  K,e  fields ,  o  is  the  numerical  grid  for  U\ 
and  A  is  the  numerical  grid  for  Ui- 


X2 


(ij-UO  s 


Fig.  4.  Control  volume ,  capital  letters  denote  the  basic  grid,  small  letters  (i,jjk)  denote  staggered  grids 
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Fig.  7.  Compact  configuration  of  four  peripheral  nozzles:  a,b,c)-a  =  5°,  d,e,f)-a  —  60°,  (a,d)-axial  velocity,  (b,e) 
concentration,  (c,f)-temperature 
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Fig,  13.  Peripheral  nozzles  are  drilled  through  the  side  walls  (Case  9).  (a,b)  axial  velocity  and  concentration  at  (x,z) 
section ,  c)  isolines  of  axial  velocity  and  vectors  (C/2,C/3)  at  cross-section  x  —  19.1ro,  d)  particle  distribution  by  temperature  at 
exit  cross-section. 
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Fig.  14.  Rectangular  cross-section  of  reactor  (Case  10).  ( a,b )  axial  velocity  and  concentration  at  (x,z)  section,  c)  isolines 
of  axial  velocity  and  vectors  (U2,Us)  at  cross-section  x  —  19.1ro,  d)  particle  distribution  by  temperature  at  exit  cross-section. 


